Collective dynamics of two-mode stochastic oscillators 
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We study a system of two-mode stochastic oscillators coupled through their collective output. 
As a function of a relevant parameter four qualitatively distinct regimes of collective behavior are 
observed. In an extended region of the parameter space the periodicity of the collective output 
is enhanced by the considered coupling. This system can be used as a new model to describe 
synchronization- like phenomena in systems of units with two or more oscillation modes. The model 
can also explain how periodic dynamics can be generated by coupling largely stochastic units. 
Similar systems could be responsible for the emergence of rhythmic behavior in complex biological 
or sociological systems. 
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Collective dynamics of stochastic systems shows a 
great variety of interesting phenomena like pulsating pat- 
terns , spiral waves |^ , synchronization |^ or stochas- 
tic resonance |^ . While the collective dynamics of single 
mode stochastic oscillators is relatively well understood 
1^ , the dynamics of stochastic oscillators with many pos- 
sible modes is less studied Q , and leads to interesting and 
new effects. An example in this sense is our recent study 
on the peculiar dynamics of the rhythmic applause [[7|, 
where the fascinating dynamics appears as a result of two 
different frequency clapping modes Many other bio- 
logical or physical systems exist, which perform stochas- 
tic oscillations with different modes. As examples we 
mention the unicellular alga Gonyaulax polyedra where 
circadian oscillators with two different periods have been 
shown to co-exist the thalamocortical relay neurons 
which can generate either spindle (7-14 Hz) or delta (0.5- 
4 Hz) oscillations and the hippocampal CAS model 
which spontaneously can generate four different rhythms 
Ipdl . When the system can switch between the available 
modes as a function of some threshold condition, new 
and interesting collective behavior is observed. 

The present study reports on such results for an en- 
semble of coupled two-mode stochastic oscillators. As 
a function of a threshold condition we reveal different 
type of synchronized and unsynchronized phases. We 
also show that rhythmic collective behavior can be gen- 
erated by suitably couphng largely stochastic units, i.e. 
the periodicity of the total output relative to the output 
of one stochastic oscillator can be strongly enhanced. As 
an immediate application we offer a new and realistic 
description for the dynamics of the rhythmic applause, 
reproducing all experimentally observed characteristics. 

Our system is composed by identical pulse-coupled 
two-mode stochastic oscillators. Their cycle can be per- 
formed in two modes: ^— i-i?/— >C— >^orA— > Bu — > 
C — > A, respectively. The periods corresponding to these 
two modes T/ and T// are given as Tj — ta + tbj + tc 



and Tji = ta + tb^ + tq , where ta , tb and tc are time 
intervals spent in state A, B and C, respectively. The 
stochastic part of the dynamics is state A, and ta is a 
stochastic variable with distribution: 

PiTA) = ^expi-^:^), (1) 

(r* — (ta))- State A should be imagined and modeled 
with an escape dynamics of a stochastic field-driven par- 
ticle from a potential valley of depth U. If the stochas- 
tic force-field is totally uncorrelated with (^) = and 
+ ''')) = DS{t) we get a distribution of escape 
times given in (1) with: r* ~ exp{U / D) . In analogy with 
the well-known FitzHugh-Nagumo system state A 
corresponds to a stochastic reaction time of the neuron 
fire. This causes all the experimentally observed fluctua- 
tions in the rhythmic human activities (l^Jl^ . In state B 
and C the dynamics is deterministic, and corresponds to 
the relaxation of the neurons [|l3|,^ . State B represents 
a "waiting time" or the rhythm giving part of the cycle. 
In biological systems this is a period the individual units 
want to impose, and usually this is the longest part of the 
cycle. The length of state B (r^^ or tbji) distinguishes 
between the two modes. We have chosen tbjj — ^tbj- 
The output of the units is in state C. During this state 
the oscillator emits a constant intensity pulse of strength 
1 /N where N is the number of oscillators in the system. 
The output of the whole system at a given moment is 

N 

f^J2f^' (2) 

1=1 

where fi is the output of oscillator i, fi = 1/N if the given 
oscillator is in state C and fi = otherwise. This total 
output is the origin of the coupling and shifts the oscil- 
lators between their operating modes. The rules for the 
evolution of the system are as follows: (i) oscillators start 
with randomly selected modes and phases and follow the 
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stated dynamics; (ii) there is a fixed output intensity, 
/*, for the system; (iii) after completing the dynamics 
in state A, each oscillator will choose to operate either 
in mode I or mode II; (iv) if at that moment / < /* 
the oscillator will operate in mode I, otherwise will fol- 
low mode II. The above dynamics has the tendency to 
keep the average total output as close as possible to /*. 
Since each oscillator has a fixed output intensity, this can 
be achieved only through switching between the available 
modes. In this sense the proposed rules are natural, mak- 
ing our model realistic. 

The evolution of a single uncoupled unit is simple since 
it preserves the starting mode (s). The average value of 
its period is: 



(T) = {ta + TB. + Tc) =T*+ TB, + TQ ■ 

Its relative standard deviation 



R = 



(T) 



T* + Tb, + TC 



(3) 
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shows an increasing tendency as a function of r* . 

In contrast with the simple uncoupled case, the dynam- 
ics of the coupled units is complex. Generally the oscil- 
lators will shift irregularly between their two modes. As 
a function of /* many interesting regimes are observed. 

For the limiting values /* = or 1 the dynamics is 
trivial. For /* = 1 we always have / < /* and all units 
operate in mode I. The / output randomly fluctuates 
due to both the initial random phases and the stochastic 
nature of ta- As the number of units increases, the vari- 
ance of this fluctuation decreases. The average value of 
the collective output (/(/* = 1)) = is a function of 
TC and (Tf) only: 



TC 

{TiY 



(5) 



The above considerations will apply in the >> 1 limit 
for the /* > cases as well. In this manner we can 
analytically predict that for high enough /* values the 
units operate as simple uncoupled stochastic oscillators 
in mode I. For /* < all units operate in mode II, and 
again there is no effective coupling between the units. 
The total signal will randomly fluctuate, and due to the 
larger (T//) value its mean will be smaller than the one 
measured in the /* > cases: 



(/>o = (/(r = 0)) = 



TC 

(Th) 



<(/)l 



(6) 



For < /* < we have a nontrivial regime, where 
the coupling is effective and the oscillators switch be- 
tween their modes. This region will be studied by com- 
puter simulations. We can numerically follow both the 
dynamics of a selected unit and the total / output of the 
system. For fixed Tg^ , r^jj and tq values the parame- 
ters governing the system dynamics are /* and r*. We 
choose tbi = 0.4 and tc = 0.1. By mapping the relevant 



{/*,T*} parameter-space, four different regimes (phases) 
can be revealed. Phase I, is an unsynchronized regime 
with (/) = (/)o (Fig. la). Phase II, is a synchronized 
regime with quasi-periodic total output and large oscilla- 
tion periods close to the (Tu) value (Fig. lb). Phase III, 
is a synchronized regime with quasi-periodic output and 
slow oscillation periods close to the (T/) value (Fig. Ic). 
Phase IV, is an unsynchronized regime with (/) = 
(Fig. Id). 

Changing the values of tb and tc leads to qualita- 
tively similar phases in the parameter space. From the 
phase-space sketched in Fig. 2 we learn that by increas- 
ing the value of /* the systems collective output changes 
from the unsynchronized phase I behavior to the synchro- 
nized phase II and phase III dynamics, leading finally to 
the unsynchronized phase IV regime. As r* is increasing 
phase III disappears (r* > 0.5) and the /* interval where 
synchronization is present (phase II and III) diminishes. 
In agreement with our previous analytic justifications, 
in phase I we have (/) — (/)o (Fig. la), in phase IV, 
(/) = and for /* > the collective behavior of 
the system always corresponds to the phase IV one. 

In the < /* < regime the dynamics of a single 
unit is non-trivial. The units stochastically shift between 
mode I and mode II oscillations. In this regime although 
the movement of each unit is largely stochastic, their col- 
lective behavior leads to a periodic output. In order to 
characterize numerically the enhancement in the period- 
icity we define a measure for it. Let us denote the output 
signal as a function of time as ,f(t). We can define an er- 
ror function, A(r), which characterizes numerically how 
strongly the f{t) signal differs from a periodic signal with 
period T 



A(T) 



2M X Jo 



f{t)-f{t + T)\dt, (7) 



where 



M = lim^^^- [ I fit)- < fit) >\ 
a; Jo 



dt 



< fit) >= lim^^oo- [ fit)dt (8) 
X Jo 

The general shape of the A(T) curve as a function of T is 
sketched on Fig. 3. For any fit) oscillating function we 
have an initially increasing tendency at small T values, 
after which for T = Tm a minimum (Am) is reached. One 
can state that Tm is the best approximation for the fit) 
signals period, and the "periodicity level" of the signal is 
characterized by: 



P = 



(9) 



We can compute this parameter both for one oscilla- 
tor working independently (pi) in the long period mode 
(where the cS'ect of stochasticity on the period is smaller) 
and for the whole system (p). The ratio p/pi will char- 
acterize the enhancement in the periodicity. We have 
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analyzed throughout the {r*,/*} parameter space this 
enhancement in the periodicity. Results for TV = 200 os- 
cillators are given in Fig. 4. From this data it is clear 
that there is a quite large parameter space (correspond- 
ing to small values of /* and large limits for r*) where the 
enhancement in the periodicity is considerable. The com- 
puted best periods (Tm) are in agreement with the phases 
previously discussed and sketched in Fig. 2. Increasing 
the number of coupled oscillators will further enhance the 
periodicity of the total output. As an example results for 
T* — 0.4 and /* = 0.1, are presented in Fig. 5. It is also 
interesting to note the stochastic resonance type effect ^] 
for the enhancement in the periodicity (Fig. 4). For fixed 
/* values there is an optimal noise level (r*) in the os- 
cillators dynamics, where a maximal p/pi enhancement 
ratio is achieved through the considered coupling. 

As an immediate application one can create a new 
model and explanation for the phenomenon of rhyth- 
mic applause [||. Different phases of the collective out- 
put (Fig. 2) reproduces the regimes observed in collec- 
tive clapping. The high intensity (phase IV) unsynchro- 
nized phase corresponds to the initial thunderous clap- 
ping of the spectators following an exceptional perfor- 
mance. Phase II describes the rhythmic applause where 
spectators clap in unison with a relatively long period. 
Phase III with a shorter period and partly synchronized 
output is characteristic for the unstable transition inter- 
vals from synchronized to unsynchronized clapping. The 
low intensity unsynchronized phase I dynamics is char- 
acteristic to the clapping of a non-enthusiastic audience. 
Modeling the clapping individuals by the proposed two- 
mode stochastic oscillators, and the coupling between 
them by the proposed threshold criteria is realistic and 
in agreement with experimental results: i) the clapping 
sound is pulse like and is reproduced by state C of the 
oscillators; ii) the clapmng period of one individual has 
a fluctuating nature |14| and is modeled by state A of 
the dynamics; iii) the rhythm of the clapping is modeled 
by state B; iv) our previously reported measurements on 
clapping individuals revealed the existence of the two 
distinct clapping modes with short and longer periods, 
respectively; v) the period of the longer clapping mode 
was found do be roughly two times larger than the one of 
the short clapping mode. In view of this new model the 
characteristic interplay between synchronized and unsyn- 
chronized regimes in the rhythmic applause should be a 
consequence of a peculiar dynamics in the /* threshold, 
and should have a psychological origin. It is clear that 
after a bad performance the /* threshold is low and leads 
to phase I type collective response. For an enthusiastic 
audience /* is big and high intensity collective response 
forms, resembling the total output of the units in phase 
IV. By lowering the level of /* (fatigue or just resting...?) 
synchronized collective response arises (phase III) which 
corresponds to the rhythmic applause. 

In our opinion the presented model has possibilities for 
understanding and modeling the origin of other rhyth- 
mic biological or sociological phenomena ll^]. The most 



interesting aspect of our results is the finding that regu- 
lar periodic output can be obtained by coupling a large 
number of stochastic units. Similar phenomena should 
be responsible for the emergence of rhythmic behavior in 
complex biological or sociological systems. 
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FIG. 1. Total / output of 200 coupled two-mode oscillators 
(r* = 0.1, TB, = 0.4, Tc = 0.1 and /* = -0.1 (a), /* = 0.1 
(b), r =0.2 (c), r=0.3 (d)). 




FIG. 2. Phases in the {/* , r*} parameter-space {N = 200, 

TBI = 0-4, TC = 0.1). 




FIG. 3. Characteristic shape of A(T). Determining the 
and Tm parameters, (r* = 0.1, /* = 0.3, tbj = 0.4, 
TC = 0.1, N = 200) 




FIG. 4. Enhancement in the periodicity (p/pi) as a func- 
tion of the T* and /* relevant parameters. (Results for 
N = 200 oscillators, tb^ = 0.4, tc = 0.1) 
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FIG. 5. Enhancement in the periodicity of the output 
(p/pi) as a function of the number of oscillators in the system. 
(/* = 0.1 and T* = 0.4 tb, = 0.4, tc = 0.1) 
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